A LAGRANGIAN APPROACH FOR THE COMPRESSIBLE 
NAVIER-STOKES EQUATIONS 



RAPHAEL DANCHIN 

■ Abstract. Here we investigate the Cauchy problem for the barotropic Navier-Stokes equations 

I in R" , in the critical Besov spaces setting. We improve recent results (see [4, 8, 9]) as regards 

the uniqueness condition: initial velocities in critical Besov spaces with (not too) negative in- 
dices generate a unique local solution. Apart from (critical) regularity, the initial density just 
has to be bounded away from and to tend to some positive constant at infinity. Density- 
dependent viscosity coefficients may be considered. Using Lagrangian coordinates is the key 
to our statements as it enables us to solve the system by means of the basic contraction map- 
ping theorem. As a consequence, conditions for uniqueness are the same as for existence, and 
Lipschitz continuity of the ffow map (in Lagrangian coordinates) is established. 
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Introduction 

We address the well-posedness issue for the barotropic compressible Navier-Stokes equations 
with variable density in the whole space M" (n > 2): 



> 

m. (0.1 



O 



dtp + div {pu) = 0, 

< dt{pu) + div {pu ®u)- 2div {n{p)D{u)) - V(A(p)divu) + V(P(p)) = 0, 

, p\\t=o = Po^ u\it=o = uo- 

Above p = p(t,x) G M+ stands for the density, u = u{t,x) £ M", for the velocity field. The 
space variable x belongs to the whole M". The notation D{u) designates the deformation tensor 



^ I which is defined by 

D{u) := ^{Du + 'Vu) with {Du)ij := djU^ and (yu)ij := diU^ . 

^ ■ The pressure function P and the viscosity coefficients A and p are given suitably smooth 

^ ! functions of the density. With no loss of generality, one may assume that P is defined over R 

and vanishes at 0. As we focus on viscous fiuids, we suppose that 

(0.2) a := min(inf (A(p) + 2p{p)), mi p{p)) > 0, 

\p>0 p>o / 

which ensures the second order operator in the velocity equation of (0.1) to be uniformly elliptic. 

We supplement System (0.1) with the condition at infinity that u tends to and p, to 
some positive constant (that may be taken equal to 1 after suitable normalization). The exact 
meaning of those boundary conditions will be given by the functional framework in which we 
shall consider the system. 

In the present paper, we aim at solving (0.1) in critical functional spaces, that is in spaces 
which have the same invariance with respect to time and space dilation as the system itself (see 
e.g. [8] for more explanations about this nowadays classical approach). In this framework, it 
has been stated [8, 9] in the constant coefficients case that, for data {po,uo) such that 

ao := (Po - 1) G B;!fm, ^0 e B;;^^^! 



Date: January 31, 2012. 



1 



2 



R. DANCHIN 



and that, for a small enough constant c, 

(0.3) II«o|IbWp(k„) < c, 

we have for any p G [1, 2n) : 

• existence of a local solution (p, n) such that a := (p — 1) G Cb{[0,T]; B^^-^), u G 

a([0,T];Bp"/r') e Li(0,r;i^;/r^); 

• uniqueness in the above space if in addition p < n. 

li p < n then the viscosity coefficients may depend (smoothly) on p and the smallness condition 
(0.3) may be replaced by the following positivity condition (see [4, 10]): 

(0.4) inf po{x) > 0. 

Those results have been somewhat extended in [16] where it has been noticed that oq may be 
taken in a larger Besov space, with another Lebesgue exponent. 

The above results are based on maximal regularity estimates in Besov spaces for the evo- 
lutionary Lame system, and on the Schauder-Tychonoff fixed point theorem. In effect, owing 
to the hyperbolicity of the density equation, there is a loss of one derivative in the stability 
estimates thus precluding the use of the contraction mapping (or Banach fixed point) theorem. 
As a consequence, with this method it is found that the conditions for uniqueness are stronger 
than those for existence. 

Following our recent paper [13] dedicated to the incompressible density-dependent Navicr- 
Stokes equation, and older works concerning the compressible Navier-Stokes equations (see [20, 
21, 22]), we here aim at solving System (0.1) in the Lagrangian coordinates. The main motivation 
is that the mass is constant along the flow hence, to some extent, only the (parabolic type) 
equation for the velocity has to be considered. After performing this change of coordinates, we 
shall see that solving (0.1) may be done by means of the Banach fixed point theorem. Hence, 
the condition for uniqueness is the same as that for the existence, and the flow map is Lipschitz 
continuous. In addition, in the case of fully nonhomogeneous fluids with variable viscosity 
coefficients, the analysis turns out to be simpler than in [4, 10] even for density-dependent 
viscosity coefficients and in the case where the density is not close to a constant. Indeed, our 
proof relies essentially on a priori estimates for a parabolic system (a suitable linearization of 
the momentum equation in Lagrangian coordinates) with rough constant depending only on the 
initial density hence time-independent. In contrast, in [4, 10] tracking the time-dependency of 
the coefficients was quite technical. 

We now come to the plan of the paper. In the next section, we introduce the compressible 
Navier-Stokes equations in Lagrangian coordinates and present our main results. Section 2 is 
devoted to the proof of our main existence and uniqueness result in the simpler case where 
the density is close to a constant and the coefficients, density independent. In Section 3, we 
treat the general fully nonhomogeneous case with nonconstant coefficients. A great deal of the 
analysis is contained in the study of the linearized momentum equation for (0.1) (see Subsection 
3.1) which turns out to be a Lame type system with variable rough coefficients. This will 
enable us to define a self-map $ on a suitably small ball of some Banach space Ep(T) and to 
apply the contraction mapping theorem so as to solve the compressible Navier-Stokes equations 
in Lagrangian coordinates. In the Appendix we prove several technical results concerning the 
Lagrangian coordinates and Besov spaces. 

Notation: Throughout, the notation C stands for a generic constant (the meaning of which 
depends on the context), and we sometimes write X < Y instead oi X < CY. For X a 
Banach space, p G [1, -|-oo] and T > 0, the notation U'{0, T; X) or Lj,{X) designates the set of 
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measurable functions f : [0,T] ^ X with t ^ \\f{t)\\x in LP(0,T). We agree that C{[0,T];X) 
denotes the set of continuous functions from [0, T] to X. 

1. Main results 

Before deriving the Lagrangian equations corresponding to (0.1), let us introduce more nota- 
tion. We agree that for a function F : ^ x R"* then divF : R'* ^ R"* with 

{divFy := diFij for 1 < j < m, 

i 

and that for A = )i<jj<n and B = (-Bjj)i<jj<„ two n x n matrices, we denote 

A : B = Tr AB = ^ AijBji. 

The notation adj(A) designates the adjugate matrix that is the transposed cofactor matrix. Of 
course if A is invertible then we have adj(A) = {det A) A~^. Finally, given some matrix A, we 
define the "twisted" deformation tensor and divergence operator (acting on vector fields z) by 
the formulae 

Da{z) := - [Dz ■ A + '^A- Vz) and div^ z ■-'^A: Vz. 

Let X be the flow associated to the vector-field u, that is the solution to 

ft 



(1.1) X{t,y) = y+ [ u{T,X{T,y))dT. 

Jo 

Denoting 



Pit,y) ■■= pit,X{t,y)) and u{t,y) = u{t, X{t,y)) 

with {p,u) a solution of (0.1), and using the chain rule and Lemma 1 from the Appendix, we 
gather that (p, u) satisfies 

f dt{Jp) = o 

(1.2) 



Jpodtu - div i^adi{DX){2n{p)DA{u) + A(p) div^ uld + P(p)Id) j = 

with J := detDX and A := {DyX)~^. Note that one may forget any reference to the initial 
Eulerian vector-field u by defining directly the "flow" X of ?2 by the formula 

(1.3) X{t,y) = y+ [ u{T,y)dT. 

Jo 

We want to solve the above system in critical homogeneous Besov spaces. Let us recall that, 
for 1 < p < oo and s < n/p, a tempered distribution u over R" belongs to the homogeneous 
Besov space S^,i(R") if 

and 

(1.4) ll"lli3;^(K") := 'Y'^^''\\'^ju\\LP(R") < oo. 

Here {Aj)j^z denotes a homogeneous dyadic resolution of unity in Fourier variables -the so- 
called Littlewood-Paley decomposition (see e.g. [1], Chap. 2 for more details on the Littlewood- 

Paley decomposition and Besov spaces). 

Loosely speaking, a function belongs to j^(R") if it has s derivatives in L^(R"'). In the 
present paper, we shall mainly use the following classical properties: 
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• the Besov space Bp i{W^) is a Banach algebra embedded in the set of continuous func- 
tions going to at infinity, whenever 1 < p < oo; 

• the usual product maps ^^(^"^M") x ^^//(R") in S^(^"^(M") whenever 1 < p < 2n; 

• Let F : / — )■ M be a smooth function (with / an open interval of R containing 0) 
vanishing at 0. Then for any s>0, l<p<oo and interval J compactly supported in 
I there exists a constant C such that 

(1-5) ^(R") - ^II^'IIb^ ^(R") 

for any a € ^(I^") with values in J. In addition, if ai and 02 are two such functions 

and s = n/p then we have 

(1-6) 11-^(02) - ^(oi)llB^/p(Kn) < C\\a2 - aiWen/p^^^y 

From now on, we shall omit M"" in the notation for Besov spaces. We shall obtain the existence 
and uniqueness of a local-in-time solution {p,u) for (1.2), with a := p — 1 in C([0, T]; 5^^^) and 
u in the space 

E^{T) := 1^; G C([0, T]; B;//"^), dtV,V^v G L^O, T; 
That space will be endowed with the norm 

Let us now state our main result. 

Theorem 1. Let 1 < p < 2n and n> 2. Let uq be a vector- field in B^^f ^. Assume that the 
initial density po satisfies oq := (po ~ 1) G ^p^i '^'^^ 
(1.7) infpo(a;)>0. 

X 

Then System (1.2) has a unique local solution {p,u) with (a, u) G C{[0,T]; B^^f) x Ep{T). 
Moreover, the flow map (ao,uo) ' — ^ {cl-iU) is Lipschitz continuous from B^^f x B^^^ ^ to 
C{[0,T];B;!f)xE,{T). 

In Eulerian coordinates, this result recasts in: 

Theorem 2. Under the above assumptions, System (0.1) has a unique local solution {p,u) 
with u G Ep{T), p hounded away from and (p — 1) G C([0, T]; B^^^). 

Let us make a few comments concerning the above assumptions. 

• We expect the Lagrangian method to improve the uniqueness conditions given in e.g. [8] 
for the full Navicr-Stokcs equations. We here consider the barotropic case for simplicity. 

• The condition 1 < p < 2n is a consequence of the product laws in Besov spaces. It 
implies that the regularity exponent for the velocity has to be greater than —1/2 (to 
be compared with —1 for the homogeneous incompressible Navier-Stokes equations). It 
would be interesting to see whether introducing a modified velocity as in B. Haspot's 
works [15, 16] allows to consider different Lebesguc exponents for the Besov spaces 
pertaining to the density and the velocity so as to go beyond p = 2n for the velocity. 

• The regularity condition over the density is stronger than that for density-dependent 
incompressible fluids (see [13]). In particular, in contrast with incompressible fluids, it is 
not clear that combining Lagrangian coordinates and critical regularity approach allows 
to consider discontinuous densities. 



COMPRESSIBLE NAVIER-STOKES EQUATIONS 



5 



• Owing to the fact that the density satisfies a transport equation, we do not expect 
Lipschitz continuity of the flow map in high norm for the Eulerian formulation to be 
true. 

• It is worth comparing our results with those of P. Germain in [14], and D. Hoff in [17] 
concerning the weak-strong uniqueness problem. In both papers, the idea is to show 
that a finite energy weak solution coincides with a strong one under some additional 
assumptions. The weak solution turns out to have less regularity as in Theorem 2. At 
the same time, the assumptions on the strong solution (p, u) are much stronger. In both 
papers, Vn has to be in -L^(0, T; L°°), and to satisfy additional conditions: roughly V^u 
or dfU have to be in L^{0,T; L'^) in Germain's work, while -/tD'^u G L'^ {0,T; L^) with 
r = 4/3 if n = 2, and r = 8/5 if n = 3 in Hoff's paper. Some regularity conditions are 
required on the density but they are, to some extent, weaker than ours. 

2. The simple case of almost homogeneous compressible fluids 

As a warm up and for the reader convenience, we here explain how local well-posedness may 
be proved for the system in Lagrangian coordinates in the simple case where: 

(1) The viscosity coefficients are constant, 

(2) The density is very close to one. 

Let fi' := A + jU. Keeping in mind the above two conditions and using the fact that the first 
equation of (1.2) implies that 

(2.1) J{t,-Mt,-) = Po, 
with J := |detDX| and 

(2.2) X{t,y) ■.= y+ [ u{T,y)dT, 

Jo 

we rewrite the equation for the Lagrangian velocity as (recall that A := {DX)~^): 

(2.3) dtu - nAu - ix'Vdivu = (1 - poJ)dtu + 2//div {&A]{DX)Da{u) - D{u)) 

+ Adiv (adj(Z)X) div^ u - divu Id) - div (adj(L»A:)P( J" Vo)) • 

The left-hand side of the above equation is the linear Lame system with constant coefficients, 
the solvability of which may be easily deduced from that of the heat equation in the whole space 
(see e.g. [1], Chap. 2). We get: 

Proposition 1. Let the viscosity coefficients (/U,/i') G satisfy > and n + p' > 0. Let 
p G [l,oo] and s G M. Let uq G B^ ^ and f G L^{0,T; B^ i). Then the Lame system 

j dtu - nAu - /x'Vdiv u = f in (0, T) x W 

\ u\t=to = Uq on 

has a unique solution u in C{[0,T)-B^i) such that dtu,V'^u G L^{0,T; B^^^) and the following 
estimate is valid: 

(2.5) ||n||^oc(B._j + min(/x,/x + /x')||V2u||^i^(^.^^) < C(||/||^i^(b|_^) + ||uo||b»^^) 

where C is an absolute constant with no dependence on fi, p! and T. 

In the rest of this section, we drop the bars on the Lagrangian velocity field. Granted with 
the above proposition, we define a map <t : v ^ u on Ep(T) where u stands for the solution to 

(2.6) dfU — juAu — /x'Vdivu = Ii{v, v) + 2iJ,div l2{v, v) + Adiv/3(v, v) — div l4{v) 
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with 

Ii{v,w) = {l-poJv)dtw, h{v,w) = ad]{DXy)DASw) - D{w), 

l3{v,w) = div A^w aAj{DX^) - divw Id, hiv) = adj{DXy)P{J~^ po). 

Note that any fixed point of $ is a solution in Ep{T) to (2.6). We claim that the existence of 
such points is a consequence of the standard Banach fixed point theorem in a suitable closed 
bah of Ep{T). 

First step: estimates for Ii, I2, I3 and I4. Throughout we assume that for a small enough 
constant c, 

(2.7) / \\Dv\\^„/pdt < c. 

Jo p.i 

In order to bound Ii{v, w), we decompose it into 

Ii{v, w) = {1 — Jy)dtw — ao(l + {Jv — ^))dtw with ao := po ~ 1- 

Taking advantage of the product law B^[\~^ x B^^f ^pli~^ (if 1 < p < 2n) and of the fact 
that Bp^f is an algebra (if 1 < p < 00), of (A.ll), (2.7) and of (1.5), we readily get 

(2.8) \\h{v,w)\\^,^^^„f,-,^ < C(||ao||^(^n/p-i) + \\dtw\\ 

Above we introduced the multiplier norm Ai{Bp^j^) for -Bp^i, that is defined by 

(2.9) II/IU(b;_j:=sup||V'/||b._^ 

where the supremum is taken over those functions tp in Bp i with norm 1. 
Next, thanks to product laws, to (A. 12), (A.13) and to (2.7), we have 

(2.10) \\l2{v,w)\\^,^^^„/,^ + ||73(«,«^)|li^(Bn/p) < Cp^;||^^(^„/,)||I)w;||^^(^W,). 

As regards the pressure term (that is /4(f)), we use the fact that under assumption (2.7), we 
have, by virtue of the composition inequality (1.5) and of flow estimates (see (A.9) and (A.ll)), 

(2-11) ll^4(?^)||^oo(Bn/P) < C(l + P?;||^^(^n/p^)(l + ||ao||^n/p). 

Second step: <1> maps a suitable closed ball in itself. At this stage, one may assert that if 
V G Ep{T) satisfies (2.7) then the right-hand side of (2.6) belongs to L^{Q,T; B'^^ ^). Hence 
Proposition 1 implies that ^{v) is well defined and maps Ep{T) to itself. However it is not clear 
that it is contractive over the whole set Ep{T). So we introduce ul the "free solution" to 

dfUL - IJ-AUL- n'VdWuL = 0, UL\t=0 = uo- 

Of course. Proposition 1 guarantees that ul belongs to Ep{T) for all T > 0. Hence, if T and 
R are small enough then any vector-field in ]3^^(j'-^[ui,R) satisfies (2.7). 

We claim that if T is small enough (a condition which will be expressed in terms of the free 
solution Ul) and if R is small enough (a condition which will depend only on the viscosity 
coefficients and on p, n and P) then 

V e Be^^t){ul,R) =^ue Be^^t){ul,R)- 

Indeed u := u — ul satisfies 

dfU — /xAS — /x'Vdivn = Ii{v, v) + 2p,d\Y l2{v, v) + Adiv/3(?;, v) — divi4(t;), 

u\t=o = 0. 
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So Proposition 1 yields^ 

Inserting inequalities (2.8), (2.10) and (2.11), we thus get: 

\\u\\e,(t) < + (ll«o||^(Bn/p-i) + \\Dv\\^,^^^„/,^)\\dtv\\^,^^^^/^^^^^ 

That is, keeping in mind that v is in Bep(t){ul, R), 

Meat) < c((l|ao||_^(^n/p-i) + P^lII^i^^^u/p^ + i?)(i? + 

So we see that, and if T satisfies 

(2.12) Cr(l + llaollon/p) < i2/2 and WDulLi tr,^/p. + WdtuiL^ (r,n/P-i. < R 

then we have 

\\u\\e,(t) < 2C(||ao||_^(^n/p-i) + 2R)R + 2CR' + R/2. 
Hence there exists a small constant r] = r]{n,p) such that if 
(2-13) ll«o||^(Bn/p-i) < r], 

and if R has been chosen small enough then u is in (j'-^{ul, R). Of course, taking R and T 
even smaller ensures that (2.7) is satisfied for all vector-field of Bej,(t){ul,R)- 

Third step: contraction properties. We claim that under Conditions (2.13) and (2.12) (with a 
smaller R if needed), the map $ is 1/2-Lipschitz over Bep[t){ul^R)- So we are given vi and 
V2 in Be^(j'^{ul,R) and denote 

ui := <l>(fi) and U2 ■= $(^2)- 

Let Xi and X2 be the flows associated to vi and i'2. Set Ai = {DXi)~^ and Jj := detDXj 
for z = 1, 2. The equation satisfied by 5u := U2 — ui reads 

dtdu — /xA&t — /x'Vdiv 5u = Sf := Sfi + 6/2 + div 5/3 + 2/xdiv (5/4 + Adiv 6/5 

with J/i := (1 - poJ2)dtdu, 6/2 := -po{J2 - Ji)dtUi, 

(5/3 := adj(L>Xi)P(poJr') - adj(L>X2)P(po Jg"'), 
(5/4 := adi{DX2)DA,{u2) - adj(DXi)L»Ai (t^i) - £'('5u), 
(5/5 := adj(DX2)^A2 : Vua - adj(r»Xi)^^i : Vui - div(5uld. 

Once again, bounding &i in Ep{T) stems from Proposition 1, which ensures that 

(2.14) ||(5u|U^(^) < (^;/-^) + r|l^/3ll.5,(^;/.) + E H^/^H^i 

i=l ' ' 1=4 

In order to bound Sfi and (5/2, we just have to use the definition of the multiplier space 
MiBplf^), and (A.11),(A.20). We get 

(2.15) < (ll«oll^(^;_/P-i) + C||Z^t;2||^^(^„/p)) ||at,5u||^^(^„/p-i^, 

(2.16) \w^\\LUB;/r') - ^ll^olU(<r')''^*''i^(</)''^*'''''i^(B;/r')- 



^For simplicity, we do not track the dependency of the coefficients with respect to and /u'. 
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Next, using the decomposition 

Sfs = (adj(DXi) - adj(DX2))P(po + adj(Z)Xi)(P(po Jf') - ^(^0^2"')), 

together with composition inequaUties (1.5), (1.6) and (A. 19), and product laws in Besov space 
yields 

(2-17) ii<5/3ii^oo(^n/P-i) <ni + ii 

Finally, we have 

,5/5 = (adj(Z?X2)-adj(Z)Xi))^^2 : Vu2+adj(DXi) ^2-^1) : Vu2 + (adj(Z)Xi)^Ai-Id) : VSu, 
whence, by virtue of (A.9), (A.IO), (A.18) and (A.19), 

(2.18) P/4|li^(^^"_/P) < ||I)*||^^(^./P)pn2||^^(^n/P) + ||Z)&.||^^(^./p^||Z)t;i||^^(^„/P). 

Bounding 6fi works exactly the same. So we see that if Conditions (2.13) and (2.12) are satisfied 
(with smaller 77 and larger C if need be) then we have 

\\Su\\Ep(T) < ^II*IIEp(T)- 

Hence, the map $ : B^^(j'-^{ul^ R) 1— t- B^^{j-^{ul,R) is 1/2-Lipschitz. Therefore, Banach' fixed 
point theorem ensures that $ admits a unique fixed point in Bep(t){ul, R)- This completes the 
proof of existence of a solution in Ep{T) for System (1.2). 

A tiny variation over the proof of the contraction properties yields uniqueness and Lipschitz 
continuity of the flow map. We eventually get: 

Theorem 3. Let p G [l,2n) (with n>\) and uq he a vector-field in B^^^ ^ . Assume that the 

initial density po satisfies ao := {po — 1) € B^^^. There exists a constant c depending only on 
p and on n such that if 

(2-19) \\^o\\m{b;{^-^) - ^ 

then System (1.2) has a unique local solution {p,u) with {a,u) G C{[0,T]; B^^f) x Ep{T). 
Moreover, the flow map (ao,«o) ' — > is Lipschitz continuous from B^^^ x B^^^ ^ to 

C{[Q,T\-B;'f)xEp{T). 

In Eulerian coordinates, this result recasts in: 

Theorem 4. Under the above assumptions, System (0.1) has a unique local solution {p,u) with 
density bounded away from vacuum and a G C([0, T]; 5^//"^) and u G Ep{T). 

We do not give here more details on how to complete the proof of Theorem 3 and its Eulerian 
counterpart. Theorem 4 as it will done in the next section under much more general assumptions. 

3. The fully nonhomogeneous case 

For treating the general case where po need not satisfy (2.19), just resorting to Proposition 1 
is not enough because the term Ii{v,v) in the r.h.s of (2.6) need not be small. One has first to 
establish a similar statement for a Lame system with nonconstant coefficients. More precisely, 
keeping in mind that p = J^^Po (we still drop the bars for notational simplicity), we recast the 
velocity equation of (1.2) in: 

Lpoiu) = h(u,u) + Pq^6xy {l2{u,u) + Iz{u,u) + h{u,u) + h{u)) 

with 

(3.1) Lp^{u) := dtu - Po ^div {2fi{po)D{u) + X{po)dWuId) 
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and 

Ii{v,w) := (1 - Jv)dtw 

h{v,w) := (adjpX,)-Id)(^(J-Vo)P^^^. + ^^.Vw;) + A(J-Vo)r^t, : Vii;)Id) 
l3{v,w) := {fi{J-^po) - n{po)){DwA^ + '^A^Vw) + (A(J-Vo) - X{po))(^A^ ■ Vw)ld 
h{v.w) := fi{po){Dw{A^, - Id) - U)Vw) + X{po)C{A^ - Id) : Vw;Id 

hiv) := -aAi{DX^)P{poJ-'). 
Therefore, in order to solve (1.2) locally, it suffices to show that the map 

(3.2) $ : t; I — >u 
with u the solution to 

J Lpoi.'^) = h{v,v) + p^'^div {l2{v,v) +l3{v,v) +I4{v,v) + l5{v)), 

[ u\t=o = uo 

has a fixed point in Ep{T) for small enough T. 

As a first step, we have to study the properties of the linear Lame operator Lp^. This is done 
in the following subsection. 

3.1. Linear parabolic systems with rough coefficients. As a warm up, we consider the 
following scalar heat equation with variable coefficients: 

(3.3) dtu - adiv (6Vn) = /. 
We assume that 

(3.4) a:= inf (a6)(t,x)>0. 

(t,a;)e[0,T]xM" 

Let us first consider the smooth case. 

Proposition 2. Assume that a and b are bounded functions satisfying (3.4) and such that bVa 
and aVb are in L°°{0,T; Bp^i) for some 1 < p < oo. There exist two constants k = k{p) and 
C = C{s,n,p) such that the solutions to (3.3) satisfy for all t £ [0,T], 

hh^^B^ ^) + '^(^M\lUb^+^) ^ (ll^ollfi.^ + j)exp('^ / ||(6Va,aV6)||^ 

whenever — mm{n/p,n/p') < s < n/p. 

Proof. We first rewrite the equation for u as follows: 

dtU — div (abVu) = f — 6Va • Vu, 

then localize the equation in the Fourier space, according to Littlewood-Paley decomposition: 

dtUj — div {abVuj) = fj — Aj(6Va • Vu) + Rj 

with Uj := AjU, fj := Ajf and Rj := div ([Aj, aftjV'u). 

Next, we multiply the above equation by Uj\uj\P~'^ and integrate over IR". Taking advantage 
of Lemma 8 in the appendix of [12] (here 1 < p < oo comes into play) and of Holder inequality, 
we get for some constant Cp depending only on p: 



n/p ' 
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which, after time integration, leads to 

(3.5) \\Uj\\i^(^LP) + CpOl2^^\Uj\\Ll(LP) < IhojIiLP 

+ \\fj\\Ll{Lp) + (l|A,-(6Va • Vu)\\lp + WRjWlp) dr. 

According to Lemmas 4 and 5 in Appendix, there exists a positive constant C and some sequence 
(cj)jez with ||c||^i(z) = 1, satisfying 

(3.6) \\Ai{bVa-Vu)\\LP + \\Rj\\Lp < Cc,-2~-^'*(||6Va|L„/p + ||aV6|L„/p) ||Vn|| . 

p,i p,i 

Then inserting (3.6) in (3.5), multiplying by 2-'* and summing up over j yields 

(3.7) IMl^^s^^^ + cpa\\u\\ ^,^^s+2^ < WuoWbs^^ + 

+ C / ||(6Va,aV6)|| r^n/p ll^ll ds+1 dr. 
Jo ^P,i 



From the interpolation inequality 
(3.8) 

we gather that 



(3-8) \\u\\j^s+i < ||u||V2 ||^i||V2 



ac 

C\\(Wa,aVb)\\.n/p\\u\\j^s+i < -^\\u\\j^s+2 + ||(6Va, aV6)|p.„/p . 

-°p,i p.i z p.i ZaCp -t>p,i p.i 

So plugging this in (3.7) and applying Gronwall lemma completes the proof of the proposition. 

□ 

In the rough case where the coefficients are only in B^^^, the above proposition has to be 
modified as follows: 

Proposition 3. Let a and b be bounded positive and satisfy (3.4) . Assume that bVa and aVb 
are in L°°{0,T; B^^f ^) with 1 < p < oo. There exist three constants rj, k and C such that if 
for some m E Z we have 

(3.9) inf(t ,r)g[o,T]x]R" Sra{ab){t, x) > a/2, 

(3.10) \\{ld-Sm){bVa,aVb)\\^^^^n/r-i^ < r?a 

then the solution to (3.3) satisfies for all t G [0,r], 

whenever 

(3.11) — mm{n/p,n/p') < s < n/p — 1. 

Proof. Given the new assumptions, it is natural to replace (3.6) by the inequality 

(3.12) \\Aj{bVa ■ Vu)\\lp + \\Rj\\Lp < Ccj2'^'{\\bVa\\-n/p-i + ||aV6|| 

which may be obtained by taking a = 1 and z/ = 1 in Lemmas 4 and 5. However, when 
bounding Rj , in addition to (3.11), one has to assume that p < n. Also, as it involves the 
highest regularity of u, we cannot expect to absorb this "remainder term" any longer, unless 
aV6 and bVa are small in B^^f~^ (which would correspond to the case that has been treated 
in the previous section). So we rather rewrite the heat equation as follows: 

dtu - div {Sm{ab)Vu) = / + div ((Id - Sm){ab)Vu) - Sm{bVa) ■ Vu - (Id - Sm){bVa) ■ Vu. 
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Now, using the infimum bound for S„i{ab) and arguing as for proving (3.5), we get 
\\uj\\L°°{Lp) +Cp"2^^lhjllLi(Lp) < Ikojllif + \\fj\\Ll{Lp) + / || Ajdiv ((Id - S'„)(a6)Vu)||LP dx 

*y 

+ ^ (||A^-(5„(6Va) • Vu)\\lp + \\Aj{(ld - Sm){bVa) ■ Vu)\\lp + \\diY {[Sm{ab), Aj]Vu)\\Lp) dr. 

The idea is to apply the procedure of the "smooth" case for the low frequency part of the 
coefficients (that is the part containing Sm) and the "perturbation" approach for the other 
part. More precisely, appealing to Lemmas 4 and 5, we get under Condition (3.11) and for some 
sequence {cj)j^z with ||c||^i(2) = 1: 

\\Ajdw {{Id -Sm){ab)Vu)\\LP < Cj2-^'\\{ld-Sm){ab)\\.r./4^u\\^s+i, 

p,l 

\\Aj{Sm{bVa) ■ Vu)\\lp < Cj2-^'\\Sm{bVa)\\^n/p\\Vu\\^s ^, 

\\Aj{{Id-Sm){bVa)-Vu)\\LP < c,-2-^--||(Id-5^)(6Va)||.„/p!i||Vt/||^.+i, 

pA 

\\diY {[Sm{ab),Aj]Vu)\\LP < Cj2-^'\\Sm'^{ab)\\-„/r,\\Vu\\^s . 

Let us plug those four inequalities in the above inequality for uj. After multiplying by 2-'* and 
summing up over j, we get 

+C(||(Id - Sm){ab)\\^^^^n/p^ + ||(Id - 5'„^)(6Va)||^^^^^„/p-i^) ||m||^i(^.+2) 

+C / \\Sjn{aVb,bVa)\\^n/p\\Vu\\r,s dr. 
Jo p>i "''^ 

It is clear that, under Condition (3.10), the second line may be absorbed by the left-hand side. 

Hence the desired inequality follows from the interpolation inequality (3.8), exactly as in the 

smooth case. □ 

We now look at the following Lame system with nonconstant coefficients: 

(3.13) dtu - 2adiv (//£>(«)) - 5V(Adivn) = /. 

Note that u and / are valued in M". We assume throughout that the following uniform ellipticity 
condition is satisfied: 

(3.14) a := min( inf {aa){t,x), inf {2aa + bX){t,x)] > 0. 

\(t,x)G[0,T]xR^ {t,x)e[0,T]xR^ J 

Let us first study the "smooth case" : 

Proposition 4. Assume that a, 6, A and [i are bounded functions satisfying (3.14) and such 
that aVfj., 6VA, jj^Va and AV6 are in L'^{0,T; B^^^) for some 1 < p < oo. There exists a 
constant C such that the solutions to (3.13) satisfy for all t G [0,T], 



)exp('-/ ||(AtVa,aV/i,AV6,6VA)||2 dr 

whenever — mm{n/p,n/p') < s < n/p. 



<C{\\uo\\bs^ + 



Proof. We introduce the following functions: 

d := |D|~^divu and := |D|^^curlM with (curlu)ij := diU^ — dju^. 

Owing to the use of homogeneous Besov space, and because the Fourier multipliers A{D) := 
|D|~^div and B{D) := |D|~-'^curl are of degree 0, it is equivalent to estimate u or {d,Cl) in 
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L'^{Bp ^) n L}p{Bp\'^). So the basic idea is to show that d and satisfy heat equations similar 
to (3.3). More precisefy, appfying A{D) to (3.13) yields 

(3.15) dtd - (2a/x + b\)Ad = A(D)(f + 2aVM • D{u) + bVX div u) 

+ [A{D),afi]Au + [A{D),afi + 6A]Vdivu. 

Given Condition (3.14), we see that arguing exactly as for proving (3.7) and because A{D) maps 
Bp I in itself, 

\\d\\L^(^Bi^)+i^a\\d\\Li^^p-2)<\\do\\ss^ \\2aV^i■ D{u) + bV\dwu\\^s^dT 

+C [ \\[A{D),a^i]Au + [A{D),an + b\]Vd\Yu\\i,s dr + C [ ||V(2a/x + 6A)|L„/p||Vu|| dr. 

Jo Jo f'l 

Note that applying Lemma 6 with a = s — 1, ly = and Lemma 4 with a = s and = yields 

||[^(D), a//]Ait|| OS < C||V(a/u)|L„/p||An|| ns-i, 

p.l p,i p.l 

\\aVn- D{u)\\^s < C||aV/i|L„/p||Vu||^, , 

p.l -°p,i p.i 

and analogous estimates for [A{D),afj, + bX]\/divu and 6VAdivu. 
Similarly, the vorticity part O of n satisfies 

dt^ - a/iAn = B{D){f + 2aV/x • D{u) + 6VAdivu) + [B{D),aiJ,]Au + [B{D), afi + 6A]Vdivn. 

So arguing exactly as for bounding d, and resorting to the interpolation inequality (3.8) and to 
Gronwall lemma, we easily get the desired inequality. It is just a matter of following the proof 
for the case of the heat equation. □ 

Let us finally focus on the "rough case" where the coefficients of (3.13) are only in L^{Bp^f). 

Proposition 5. Let a, b, A and /x be bounded functions satisfying (3.14). Assume that dV/J,, 
6VA, /uVa and XVb are in L°°(O..T; B^^f ^) for some I < p < oo. There exist two constants 
r] and n such that if for some m G Z, we have 

(3.16) min(^inf(t^^)g[o,r]xMn Smi^an + 6A)(t, x), inf(t^^)g[o,r]xM" Smian){t,xjj > 

(3.17) ||(Id- S^)(/iVa,aVM,AV6,6VA)||^^^^„/P-i) < va 
then the solutions to (3.13) satisfy for all t £ [0, T], 

\MLoo(B;^^)+(^\MLj{B;f) 

<C(||no||^. +11/11^1(5. ))exp(- / ||5^(/iVa,aV/x,AV6,6VA)f /pdr) 
p.l ty p,i) \a Jq tSp,i J 

whenever (3.11) is satisfied. 

Proof. As for the heat equation, wc split the coefficients of the system into a smooth (but 
large) low frequency part and a rough (but small) high frequency part. It turns out to be more 
convenient to work directly on the equations for d and VL. More precisely, as regards d, we write 
(starting from (3.15) and denoting c := 2a/x + 6A) that 

dtd - div (cVd) = -Vc • Vd + A{D){f + 2aV ii ■ D{u) + bVXdwu) 

+[A{D),an]Au+ [^(Z)),a/x + 6A]Vdivn, 



a 
2' 
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whence, denoting dj := Ajd, 
dfdj — div (SmcVdj) = div {[Aj, Smc\Vd) 

+Aj (div ((Id - Srn)cVd) - SmVc ■ Vd - (Id - S'„)Vc • Vdj 
+AjA{D){f + 2Sm{aVn) ■ D{u) + 2(Id - Sm){aVii) ■ D{u) 
+Sm{bVX) dwu + {Id- SmKbVX) div u) + Aj([A{D), Sm{an)]Au 

+ [A{D), Smian + 6A)]Vdivn + [A{D), (Id - Sm){aiJ,)]Au + [A{D), (Id - 5^)(a/x + 6A)]Vdivn) . 

Under Condition (3.11), Lemmas 4, 5 and 6 imply that 



||div([A^-,5„c]Vd)||LP 


< 


c,2- 


-^■^||5™Vc||^„/p||Vd||^.^, 
-^■^||(Id-5jjc||^„/,||Vd||^.+i, 


||Ajdiv((Id - Sm)cVd)\\LP 


< 


c,2- 


\\Aj{SmVc-Vd)\\LP 


< 


c,2- 


-^■^||5™Vc||^„/,||vi||^.^, 

-^■^Kld - ^J)Vc|U„/p-i||Vd||^.+i, 


\\Aj{{ld- Sm)Vc-Vd) 


< 


c,2' 


\\Aj[A{D),S^{afi)]Au\\LP 


< 


c,2- 


-^■''\\VSm{afJ,)\\ ^n/p \\Au\\ , 
-i^||V(Id - .„/,!i||Atx||^. 


\\Aj[A{D), {Id - Sm){afi)]Au\\LP 


< 


Cj2- 



and similar estimates for 

AjA{D){Sm{aVfi) ■ D{u)), A,.4(L>)((Id - -S^)(aV/x) • D{u)), 
AjA{D){Sm{bVX) div u), AjA{D){{ld - Sm){bVX) div u), 
Aj[A{D),Sm{an + bX)]Vdivu, Aj[A{D), (Id - Sm){an + bX)]Vdivu. 
Of course the curl part Q of the velocity may be treated in the same way. Therefore we get 



+ / ||<S'm(aV//,//Va, 6VA, AV6)|L„/p||'u|| ns+i dr 

Jo P'l ^'^ 

+ ||(Id-5„)(aV/x,MVa,6VA,AV6)|| . dr. 

Jo P''^ 



Obviously the last term may be absorbed by the left-hand side if r] is small enough in (3.17) and 
the last-but-one term may be handled by interpolation according to (3.8). So applying Gronwall 
lemma yields the desired inequality. □ 

For the sake of completeness, we still have to justify the existence of a solution to (3.13). 
More precisely, we want to establish the following result: 

Proposition 6. Let p be in (l,+oo). Let a, b, A and fi be bounded functions satisfying (3.14). 
Assume in addition that there exist some constants d, b, A and p, such that 

(3.18) 2afl + bX>0 and dfl > 0, 

and such that a — d, b — b, fJ, — and A — A are in C([0, T];Bp^f). Finally, suppose that 

(3.19) lim \\{Id-Sm){a-d,b-b,X-X,n-j2)\\^^„/p=0. 

Then for any data uo G B^ ^ and f G L^(0, T; j^) with s satisfying (3.11), System (3.13) 
admits a unique solution -u G C([0, T]; i^^ ]^)nL^(0, T; Bp'^). That solution satisfies the estimates 
of Proposition 5 for all large enough m E 7,. are fulfilled. 
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Proof. The proof is based on the continuity method as explained in e.g. [18] (and used in [11] 
in a similar context as ours). For 6 G [0,1], we introduce the following second order operator 
Pg acting on vector-fields u as follows: 

Vqu := — 2a0div {iX0D{u)) — 6gV(Aediv'u), 

where (iQ := (1 — 9)a + 9a, bg := (1 — 9)b + 9b, and so on. We claim that one may find some 
m £ Z independent of such that for all 9 G [0, 1], the conditions (3.16) and (3.17) are satisfied 
by ag, bg, fig and \g. Indeed, we notice that 

ag — a = 9{a — a). 

Hence, for all 9 e [0, 1], 

||(Id - Sm){ag - a)ll^oc(^"/P) < ll(Id - Sm){a - 

and similar properties hold for bg, Xg and /ig. In particular, owing to the continuous embedding 
of Bp^f in the set of continuous bounded functions, and to (3.19), wc deduce that there exists 
some m e Z so that the elhpticity condition (3.16) is satisfied by operator Vg for all 9 G [0, 1]. 
Likewise, we have for instance 

^igVag = 9{1 - e)fiVa + e'^/j.Va 

and similar relations for the other coefficients. Hence one may find some large enough m so 
that (3.17) is satisfied for all 9 G [0, 1]. In addition, the above relation shows that 

||S'm(Ai9Vae)|Ln/p < /i||S'^Va||^n/p + \\Sm{l^Va)\\^n/p. 

p,i p,i p,i 

Hence all the terms appearing in the exponential term of the estimate in Proposition 5 may 

be bounded by a constant depending only on m and on the coefficients a, b, A and fi. As a 
conclusion, one may thus find some constant C independent of such that any solution w of 

dtw - Vgw = g, w\t=o = wq 

satisfies 

(3.20) \H\L^iBl,)+<^\H\LUB;f) ^ C-dl^ollB.^^ + MIl^B^^ ,))- 

After this preliminary work, one may start with the proof of existence (uniqueness follows from 
the estimates of Proposition 5). Let £ be the set of those 9 in [0,1] such that for every data 
uq and / (as in the statement of the theorem), System 

(3.21) dtu - Vgu = f, u\t=o = uo 

has a solution u in the set F^{T) := C([0, T]; J n L^{0, T; B^'^i^). 

Note that according to Proposition 1, the set £ contains hence is nonempty. So it suffices 
to find a fixed e > such that for all 9o e 6, we have 

(3.22) [9o-e,eo + e]n[0,l]cS. 

So let us fix some ^ £, uo ^ B^ i, f G ^''^(0, T; B^ i) and v G F^iT) and consider the solution 
u to the system 

dtU - Vg^U = f + {Vg - Vgo)v 
with 9 G [0,1] such that \0 — 9o\ < e. Given that is in £, the existence of u in Fp(T) is 
granted if {Vg — Vgg)v G L^{0, T; B^^^). So let us first check this: we have 

{Vg - Vgo)v = {9- ^o) (2aeodiv {{fi - n)D{v)) + 2{a - a)div {ngD{v)) 

+6^0 V ((A - X)dWv) + (b- b)V{XgdWv)y 
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Under Condition (3.11), one may thus conclude thanks to product estimates in Besov spaces 
(see Lemma 4) that {V0 - Ve^jv G L^(0,r; 5^ 1). Furthermore 

\\iVe - ^60)^^115" < Ce[{a + Wae^ - a|L„/p)||^ - /i|Ln/p + (A + WlJ^e - AllR"/p)l|a - a|Ln/p 

+(6+ 116^0 - ^IIr"/p)I|A - ^llon/P + (A+ \\\e - A|| •„/p)||a - a|| -n/p J 

The coefficients may be bounded in terms of the initial coefficients a, 6, A and ii. Hence, 
applying (3.20) we get for some constant independent of Oq and of 0, 

Taking e small enough, it becomes clear that the linear map : v ^ u \s contractive on the 
Banach space F^{T). Hence it has a (unique) fixed point u G F^{T). In other words, u satisfies 
(3.21). 

Given that £ is nonempty and that e is independent of one may now conclude that 1 is 
in £. Therefore, there exists a solution u G F^{T) to (3.13). □ 

Remark 1. Under the assumptions of the above proposition, the constructed solution u satisfies 
dtu G L^{0,T; Bp^^). Indeed, it suffices to notice that 

dfu = f + {a + {a- a))div {p. + {fi - p)D{u)) + (6 + (6 - 6))V(A + (A - A)divu), 

and to use Lemma 4 together with the facts that Vu is in L^(0,r;5*+^). Moreover we have 

ll^t«llLi(B^) <C'(||«o|Ib.+ 11/11 rVB.Oexpf^ / \\Sm{^iVa, aV 1^, XVb,bV X)f dr) 

J PA' PA t ^ PA' \ Ct JQ P,l / 

where C may depend also on the norm of a — a, b — b, A — A and fi — jl in L°°(0, T; B^^^). 

3.2. Proof of Theorem 1. As we want to consider (possibly) large velocities, we introduce, as 

in the almost homogeneous case the free solution to the Lame system corresponding to p = 1, 
that is the vector-field ul in Ep{T), given by Proposition 1), satisfying^ 

LiUL = 0, u\t=Q = Uq. 

We claim that the Banach fixed point theorem applies to the map $ defined in (3.2) in some 
closed ball BE^^rp){uL,R) with suitably small T and R. Denoting u := u — ul, we see that u 
has to satisfy 

I LpoU = Ii{v,v) + pQ^div {l2{v,v) + l3{v,v) +I4{v,v) +/5(^^)) + (-^1 - Lpg)uL, 
I u\t=o = 0. 



(3.23) 



If the right-hand side is in L^{0,T; B^^^ ^) and if there exists some m G Z so that (3.16) and 
(3.17) are fulfilled then Proposition 6 and Remark 1 ensure the existence of u in Ep{T). Now, 
the existence of m so that 



mm 



( inf 5^(2^ + ^\msJ^\] > " 



and 



I — iiv 1 — ' 17 iiv I I I r\ 

\x \ po po^^ ^ Po ^ J 2 

(Id - S^) (^Vpo, ^Vpo, ^Vpo, ^Vpo) 
^ Po Po Po Po ^ 



. , , < ria. 

An/p-l — ' 



is ensured by the fact that all the coefficients (minus some suitable constant) belong to the space 
Bp^f the definition of which is gi 
(3.23) will be carried out below. 



Bp^f the definition of which is given by a convergent series. The study of the right-hand side of 



See (3.1) for the definition of operator L. 
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First step: Stability of Be^(j^{ul, R) for small enough R and T . Proposition 5 and the defini- 
tion of the multipher space M{Bp^f~^) ensure that 

(3.24) \\u\\e,{t) < Ce^^O'"'^(||/i(t;,i;)||^i^(^n/p-i) + - Lp,)uL\\^,^^^n/p-i^ 

for some constant Cpa^m depending only on pQ and on m. 

In what follows, we assume that T and R have been chosen so that (2.7) is satisfied by v. 
From (A. 11) and product estimates, we have Iiiv,w) £ L^iO,T; B'^^f ^) if also w is in Ep{T), 
and 

(3.25) \\h{v,w)\\^,^^^„j,-^^ < 

Next, using the decomposition 

(Li - Lpq)ul = (Po ^ - l)div {2ii{pq)D{ul) + \{pQ)divuL Id) 

+div (2(^(po) - l^{l))D{u) + (A(po) - A(l))div«Id), 

and composition inequalities (1.5) and (1.6), we see that (Li - Lp,,)uL G L^{0,T;Bplf ^) and 

(3.26) ||(Li - LpJ«L||^i^^^„/j,_i^ < ||ao||^n/p(l + ||oo||^n/p)||i:'uL||^i^^^„/p^. 

Likewise, flow and composition estimates ensure (under assumption (2.7)) that 

(3.27) \\Ii{v,w)\\^,^^^n/r,-i^ < (1 + \\ao\\^n/r)\\Dv\\^,^^^^/,^\\Dw\\^,^^^^^ for z = 2,3,4 

and that 

(3.28) ll-^5(^)||^^(3n/p) < r(l + llaoll ^„/p)(l + \\Dv\\^,^^^^/,^). 

So plugging the above inequalities in (3.24) and keeping in mind that v to satisfies (2.7), we get 
after decomposing u into u + ul: 



Meat) < Ce^""-™' (1 + ||ao||^n/p)' (^(T + ||ao||^n/p||I?nL||^i^(^„/p^) 
So, because v £ B^^^^j^-^^ul, R), 

Meat) < Ce^^O'^^il + ||ao||^„/p)'((T + ||ao||^„/p||L»tXL||^i^(^„/pP 

Therefore, if we first choose R so that for a small enough constant rj, 

(3.29) (l + ||ao|Lwp)'i?<r? 

and then take T so that 

Cpo,mT < log 2, T < llaoll WpP^^lILi .^n/p. < 

(3.30) '"^ ^ 

then we may conclude that $ maps Bej,(t){ul, R) into itself. 
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Second step: contraction estimates. Let us now establish that, under Condition (3.30) $ is 
contractive. We consider two vector-fields and in Bej,[t){'^l, R) and set := ^{v^) and 
:= Let Su := — and 5v := v'^ — v^. In order to prove that $ is contractive, it is 

mainly a matter of applying Proposition 5 to 

Lp> = h{v\Sv) + (J^i - J„2)att;2 + p-Miv {{hiv^v^) - l2{v\v^)) 

+{h{v\v^) - h{v\v')) + (hiv^v^) - h{v\v')) + ihiv^) - h{v')))- 
So we have, given that Cp^^mT < log 2, 

(3.31) \\5a\\ Ep{T) < C(l + l|ao||^n/p)(||/i('y\&;)||^i^^^„/p-i^ + ||(J^i - J^2)dtv'^\\^,^^^r./p-i^ 

The first term of the right-hand side may be bounded by means of (3.25). As for the second 
term, product estimates and (A.20) imply that 

||(J„i - J^2)dtv'^\\^i^^^n/p~i^ < C\\Ddu\\^i^^^„/p^\\dtv'^\\^^^^„/p-iy 

In order to deal with the next term, we use the decomposition 

hiv^v^) - h{v\v^) = A(J-2Vo)r^2 : Vt;2)(adj(DX„2) - adj(DX,i)) 

+ (adj (15X^1 ) - Id) (A(J-2Vo) - A(J-iVo)) (^^2 : V^;^) 

+ (adj(DX,i) - Id)A(J-,Vo)((^^2 - ^A„i) : Vv^ + : V5c) 

+ terms pertaining to /x. 

Taking advantage of product laws in Besov spaces, of composition estimates (1.5) and (1.6), and 
of the flow estimates in the appendix, we deduce that for some constant Cp^ depending only on 

\\h{v\v^)-h{v\v^)\\^,^^^„jp^<Cp,\\{Dv\Dv^^^^ 

Similar estimates may be proved for the next two terms of the right-hand side of (3.31). Con- 
cerning the last one, wc use the decomposition 

h{v^) - h{v^) = (adj(i)X,i) - adj(L>X,2))P(J-,Vo) - adj(L>X,i)(P(J-,Vo) - ^'(^-iVo))- 
Hence 

\\h{v^) - h{v^)\\^,^^^;lr^ < C{1 + ||ao||^„/.)r||i?<5u||^^^^„/,^. 

We end up with 

\Mep(t) < C(l + ||ao||^„/p)'((r+ ||(L't;\£'^;')||^i^(^n/p) + ||at^;2||^^^^„/,_ip||Z)&||^^^^„/,^ 

Given that and v"^ are in Be^(j,){ul,R)-, our hypotheses over T and R (with smaller rj in 
(3.29)) thus ensure that, say, 

WHEpiT) < \\\Hep{T)- 

One can thus conclude that $ admits a unique fixed point in B^^(^j^-^{ul, R). 
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Third step: Regularity of the density. Granted with the above velocity field u in Ep(T), we set 
p := Ju^po- By construction, the couple (p, u) satisfies (1.2). Let us now prove that a := p — 1 
is in C([0,r];Bp_(^). We have 

Given (A. 11) and using the fact that Du G L^{0,T; B^^f), it is clear that Ju^"^ belongs to 

C{[0,T]; Bp^f). Hence a belongs to C{[0,T]; B^^f), too. Because B^^^ is continuously embedded 
in L°°, Condition (1.7) is fulfilled on [0,T] (taking T smaller if needed). 

Last step: Uniqueness and continuity of the flow map. We now consider two couples {pq,Uq) 
and (pqjUq) of data fulfilling the assumptions of Theorem 1 and we denote by {p^,u^) and 
(p'^jU^) two solutions in Ep{T) corresponding to those data. Setting du := — u^, we see that 

Lpli^) = {^pI - -^p2)(^i^) + (1 - Ju^)dtSu + (J„i - Ju2)dtu^ 

+{plr'd\y (^((/| - I}){u\u') + (/| - Il){u')), 

where Ig, /|, II and /g correspond to the quantities that have been defined just above (3.2), 
with density Pq. Note that those terms and also (1 — Jui)dt5u and (J„i — Jy2)dtu^ may be 
bounded exactly as in the second step. So the only definitely new terms are (L^i — Lp2){u'^) 

and the last line. As regards (L^i — L^2)(n^), it may be decomposed into 

{Lp.-L^){u'') = ((pi)-i-(p2)-i)div(2^(pi)i^(n2) + A(pi)div^2id) 

- (pg)-Miv {2{M) - M))D{u') + WpD - X{pl))divuHd). 
Hence, combining composition, flow and product estimates, we get for t <T, 

It is not difficult to show that the other "new" terms satisfy analogous estimates. Hence, 
applying Proposition 5 to the system that is satisfied by 6u, we discover that for t <T, 

+ + \\^0\\^n/p + \\Spo\\^n/p{t+ ||(W, W)||^i^^„/p^)). 

Let us emphasize that the constant Cpi,p2 depends only on Pq through its norm for the integer 
m used in Proposition 5 corresponds to pg only. Hence if 5po is small enough then the above 
inequality recasts in 

+ ll&to|Ln/p + \\Spo\\^u/p{t + \\DU^\\^ -^/p )). 

An obvious bootstrap argument thus shows that if t, 5uq and 5po are small enough then 

\\5a\\E^(t) < 2Cpo(ll*^o||o"/F + WoWr^n/p). 
As regards the density, we notice that 

da = J-,^Sao + {J-^ - J-^)al. 
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Hence for all t G [0,r], 

WHMsn/p < C(l + \\DU^\\^ -n/pM&loW^n/pWDSuW^^n/p . 

p,l t \ p,l ' p,l t V p,l / 

So we eventually get uniqueness and continuity of the flow map on a small enough time interval. 
Then iterating the proof yields uniqueness on the initial time interval [0, T]. Note that it also 
yields Lipschitz continuity of the flow map for the velocity as for fixed data {pq,Uq), one may 
find some neighborhood and common time interval on which all the solutions constructed in the 
previous steps exist. 

3.3. Proof of Theorem 2. For uq G B^^f~^ and po £ (1 + -Sp^i^), the local existence for 

(0.1) may be proved directly (see [4, 8]) hut only under the assumption that p < n in the case of 
nonconstant viscosity coefficients. Here we get the result (including uniqueness) from Theorem 1, 
and under the sole assumption that p < 2n. This is a mere corollary of the following proposition 
which states the equivalence of the systems (0.1) and (1.2) in our functional setting. 

Proposition 7. Assume that the couple {p,u) with {p — I) G C{[0,T]; B^^f) and u G Ep{T) 
(with 1 < p < 2n) is a solution to (0.1) such that 

(3.32) / \\Vu\\-n/pdt<c. 

Jo ^p.i 

Let X he the flow of u defined in (1.1). Then the couple {p,u) := {poX,uoX) helongs to the 
same functional space as {p,u), and satisfies (1.2). 

Conversely, if {p— l,u) helongs to C{[0,T]; B^^f) x Ep{T) and {p,u) satisfies (1.2) and, for 
a small enough constant c, 

(3.33) / \\S/u\\^n/pdt<c 

Jo f.i 

then the map X defined in (1.3) is a (and in fact a locally BjJi~^^) diffeomorphism over M."' 
and the couple {p,u) := {po X~^ ,uo X~^) satisfies (0.1) and has the same regularity as {p,u). 

Proof. Let us first consider a solution {p,u) to (0.1) with the above properties. Then, the 
definition of X implies that DX — Id is in C([0, T]; i?p(^). In addition. Proposition 8 ensures 
that {p,u) := (p o X,u o X) belongs to the same functional space as {p,u), and (A. 9), (A. 10), 
(A.ll) below imply that A - Id, adj(L>X) - Id and J'^ - 1 are in C([0, T]; 5p(^). Therefore 
the product laws for Besov spaces enable us to use the algebraic relations (A. 5), (A. 6), (A.7) 
and (A. 8) whenever p < 2n. Therefore {p,u) fulfills (1.2). 

Conversely, if we are given some solution (p, u) in C([0,T]; (1 + Bplf)) X Ep{T) to (1.2) then 
one may check (see the appendix of [13]) that, under condition (2.7), the "flow" X{t,-) of u 
defined by 

(3.34) X{t,y):=y+ fv{T,y)dT 

Jo 

is a diffeomorphism over M", and satisfies DX - Id G Ci[Q,T]; B^jf). Hence one may 
construct the Eulerian vector-field u and Eulerian density by setting 

p{t,-) := po X~'^{t,-) and u{t,-) ■= uo X~'^{t,-)- 

As above, the algebraic relations (A. 5), (A. 6), (A.7) and (A. 8) whenever p < 2n, hence {p,u) 
is a solution to (0.1). That {p,u) has the desired regularity stems from Proposition 8. □ 
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Proof of Theorem 2. We consider data (po^^^o) with po bounded away from 0, (po — 1) G 
Bp^f and uo G ^ Then Theorem 1 provides a local solution {p,u) to System (1.2) in 

C([0,r];(l + Bplf)) X Ep{T). If T is small enough then (3.33) is satisfied so Proposition 7 
ensures that (p o X~^,u o X"^) is a solution of (0.1) in the desired functional space. 

In order to prove uniqueness, we consider two solutions {p^,u^) and {p^^v?) corresponding 
to the same data {pq,uq), and perform the Lagrangian change of variable (pertaining to the 
flow of li^ and v? respectively). The obtained vector-fields and v? are in Ep{T') and both 
satisfy (1.2) with the same po and uo- Hence they coincide, as a consequence of the uniqueness 
part of Theorem 1. □ 



Appendix A. Appendix 

A.l. Change of coordinates. Here we establish a result of regularity concerning changes of 
variables in Besov spaces. Even though this kind of result is somewhat classical (at least in 
nonhomogeneous Besov spaces), we did not find any reference in the literature of the estimates 
that we need. We here give a result in general Besov spaces Bp ^, the definition of which may 
be found in e.g. [1] . 

Proposition 8. Let X be a globally bi-Lipschitz diffeomorphism of and {s,p,q) with 1 < 
p < oo and —n/p' < s < n/p (or just —n/p' < s < n/p if q = 1 and just —n/p' < s < n/p if 
q = oo). 

Then a i-^ ao X is a self-map over Bpg in the following cases: 

(1) SG{0,1), 

(2) s G (—1,0] and Jx-i ^he multiplier space •M.{B^,^^,) defined in (2.9), 

(3) s > 1 and {DX - Id) G 

Proof: Let us first assume that s G (0, 1) and q = p. Then one may use the classical charac- 
terization of the norm of Bp p in terms of finite differences (see e.g. [1]) so as to write: 

ff f HX(y))-u(X(x))\P , , 
; \Jf.nJf,n \y - 

Hence performing the change of variable x' = X{x) and y' = X{y), we get 



1 



whence 

The condition that s < n/p ensures in addition that u belongs to some Lebesg ue space Lp* (M ) 
with p* < oo (or in the set of continuous functions going to at infinity if g = 1 and s = n/p). 
Hence uoX G Lp*(]R") too and one may thus conclude that uoX ^ Bp^p{W^). An interpolation 
argument then yields the desired result for any s G (0, 1) and q G [l,-\-infty]. 

The result for negative s may be achieved by duality: we have 

||moX||o<, (mn)= sup / v{z)u{X{z))dz. 

"'"^ ' \\v\\ -s <iJm.^ 

II lis / ,(H")- 
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Now, setting x = X{z), we have 

/ v{z)u{X{z))dz = [ u{x)v{X-^{x))dx, 

= / u{x)v{X''^(x))Jx-i{x) dx, 

^ II^IIb|,,(R")II^ ° X'^Jx-^ IIb;/,,(m")- 

So the definition of the multiplier space and the first part of the lemma allows to conclude. 

Finally, let us examine the cases of larger values of s. If 1 < s < 2 then one may write 

D{uoX) = (DuoX) -DX. 

AsO<s — 1<1, the first part of the proof ensures that Du o X G ^p^q^- moreover 

{DX — Id) G B^^i , the standard product laws in Besov spaces give the result. 
If 2 < s < 3 then we use the algebraic relation, 

D'^{u oX) = {D^u o X){DX, DX) + D^X ■ {Du o X). 

Hence the result follows from product laws and the previous result applied with s — 1 or s — 2. 

The higher values of s may be achieved by induction, and the remaining cases (s an integer) 
follow by interpolation. The details are left to the reader. ■ 

A. 2. Some properties of Lagrangian coordinates. Let us first derive a few algebraic re- 
lations involving changes of coordinates. We are given a -diffeomorphism X over M". For 
H : ^ W^, we agree that S{y) = H{x) with x = X{y). With this convention, the chain 
rule writes 

(A.l) DyH{y) = D,H{X{y))-DyX{y) with {D^H)ij = d^.W and {DyX)ij = dy.X\ 
or, denoting Vy = ^Dy and V^; = ^Dx, 

VyH{y) = {VyX{y)) • V,H{X{y)). 

Hence we have 

(A.2) DxH{x) = DyH{y)-A{y) with A{y) = {DyX{y))-' = D^X-^x). 

Lemma 1. Let K be a scalar function over M" and H, a vector-field. Let X be a 
diff'eomorphism such that J := det{DyX) > 0. Then the following relations hold true: 



(A.3) Va^K = dWy {adi{DyX)K), 

(A.4) divxH = J-^ divy {adi{DyX)H), 

where ad}{DyX) stands for the adjugate of DyX. 

Proof: The first item stems from the following series of computations (based on integrations 
by parts, changes of variable and (A.2)) which hold for any vector-field (p with coefficients in 
C~(M"): 

jv^K{x).mdx = - j K{x)d-^x)dx, 

= -J K{y)divx(f){y)J{y)dy 

= -J J{y)K{y)Dy4>{y):A{y)dy, 

= J 4>{y) ■ divy {adj{DyX)K){y) dy, 

= j ^{x) - diYy{ajd]{DyX)K){X-^{x))J~^{X-^{x))dx. 
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Proving the second item is similar. ■ 
Combining (A.2), (A.4) and (A.3), we deduce that if u : M" ^ M" and P : M" ^ R then 

(A.5) ~S~a = dWy (adj {DyX)V~^) = diVy (adj {DyX)'^AVyu) , 

(A.6) Vx divj; u = J"^ dWy {adi{DyX) div^ u) = J"^ divy {ajdj{DyX)'^A : Vyu), 

(A.7) = divj, {adi{DyX)P). 

Note that we will use the above relations in the case where X is the flow of some time-dependent 
vector field u, defined by the relation 

Xit,y)=y+ [ u{T,X{T,y))dT for all te[0,T]. 
Jo 

Hence we will also Imvc 



(A. 8) dtp + div (pu) = dt{Jp) and dt{pu) + div {pu ®u) = dt{Jpu). 

Let us now establish some estimates for the flow Xy of some given "Lagrangian" vector field 
(that is Xy is defined by (3.34)). 

Lemma 2. Let p G [1,+oc) and v be in Ep{T) satisfying (2.7). Let Xy be defined by (3.34). 
Then we have for all t G [0,T], 

(A.9) ||Id - adj(L'X„(i))|U„/. < \\Dv\\ 

(A.IO) ||Id - Am^n/p < \\Dv\\ , 

(A.ll) \\JtHt) - 111 < \\Dv\\^,^^„/,y 

Furthermore, if w is a vector field such that Dw G L^(0,T; B^^^) then 

(A.12) \\{adi{DXy)DAAw) - Diw))it)\\^^/, < \\D^^,^^^/,^\\Dw\\^,^^^/py 

(A.13) ||(adj(r>X^)divA,(u;) - divu) Id)(t)||^„/p < ||r>i;||^i^^„/p^||I?w;||^i^^„/p^. 

Proof: Recall that (see e.g. the appendix of [13]) for any n x n matrix C we have 
(A. 14) Id - adj (Id + C) = {C - (Tr C7)Id) + P2{C), 

where the entries of the matrix P-ziC) are at least quadratic polynomials of degree n — 1. 
Applying this relation to the matrix DX(t), and using the fact that 

(A.15) DXy{t,y)-ld= ['Dv{T,y)dT, 

Jo 

we deduce that 

Id - adj{DXy{t)) = j^[Dv- div Id) dr + P2 (^(^ Dv dr)^ . 

Given that is a Banach algebra and that (2.7) holds, we readily get (A.9). 

In order to prove (A.IO), we just use the fact that, under assumption (2.7), we have 

(A.16) Ay{t) = {ld + Cy{t))-^ = Y,{-'^f{Cy{t)f with Cy{t)= t DvdT, 



and that is a Banach algebra. 
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As regards (A. 11), we write 

(A.17) Mt,y) = l+ [ dwv{T,X^{T,y))J^{T,y)dT = l+ [ {Dv : Bdi{DX^)){T,y) dr. 

Jo Jo 

Hence, if Condition (2.7) holds then we have (A. 11) for J. In order to get the inequahty for 
J^^, it suffices to use the fact that 

J-\t,y) - 1 = (1 + {Mt,y) - 1))-^ - 1 = V(-l)'= fov : adj(L>X^)cir. 

k>i 

For proving (A. 12), we use the decomposition 

2{sidi{DX^)DAAw) - D{iB)) = (adiiDX^) -ld)iDw + Viu) 

+(Id + (adj(DX^) - Id)) {Dw ■ {A - Id) + (^A - Id) • Vw) . 

Hence the desired inequality stems from (A. 9) and (A. 10), and from the fact that B^^f is a 
Banach algebra. Inequality (A. 13) is similar. This completes the proof of the lemma. ■ 

Lemma 3. Let vi and V2 be two vector-fields satisfying (2.7), and Sv := V2 — vi. Then we have 
for all p £ [l,+oo) and all t G [0, T] (with obvious notation): 

(A.18) \\A2{t)-Ai{t)\\-nj, < p&ll , .„/p., 

(A.19) ||adj(DX2(i)) - adj(i^Xi(i))||^n/. < ||L>&||^,^„/, 

(A.20) \\4\t) - Jt\t)\\ ^„/p < \\Ddv\\^,^^^^/,y 

Proof: In order to prove the first inequality, we use the fact that, for i = 1, 2, we have 

rt 



Ai = (Id + Ci)~^ = V(-l)'=Cf with Ci{t) = f Dvi dr. 
k>o 



Hence 

a2-a,=y: (ci - ) = (fDdu dr) y: e CiCt'-'. 

k>l ^-^^ ^ k>lj=0 

So using the fact that Bp^f is a Banach algebra, it is easy to conclude to (A.18). 

The second inequality is a consequence of the decomposition (A. 14) and of the Taylor formula 
which ensures that, denoting SC := C2 — C\, 

adj(DX2) - adj(DXi) = (Tr((5C))Id - dC + dP2{Ci){dO) + ^d^ P2{Ci,Ci){dC, dO) + ■ ■ ■ 

where the coefficients of P2 are polynomials of degree n — 1. As the sum is finite and is a 

Banach algebra, we get (A.19). 

Proving the third inequality relies on similar arguments. It is only a matter of using (A.17). 
The details are left to the reader. 
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A. 3. Commutator and product estimates. This last paragraph is devoted to the proof of 

commutator and product estimates that have been used for investigating the Lame system. 
Those proofs rely on the following Bony decomposition (first introduced in [2]) for the product 
of two functions: 

(A.21) fg = Tfg + R{f,g)+Tgf. 

The paraproduct and remainder operators T and R are defined by 

Tfg:= ^ Aj,fAjg and R{f,g) := Yl ^j'f^id, 
j'<j-^ \j'-j\<i 
where (Aj)^^^ stands for some homogeneous Littlewood-Paley decomposition. 
Lemma 4. Let p G [1, +oo] and the real numbers v and a satisfy 



1/ > and — min 



f n n\ n 

-,- <<T< 
\P p'J p 



Then the following estimate holds true for all tempered distributions f and g over MP : 

Proof. The result relies on Bony decomposition (A.21). The standard continuity results for the 
paraproduct and remainder operators ensure that (see e.g. [1], Chap. 2): 

\\T{g\\B- ^\\f\\B--\\g\\B''+- ifz/>o, 

00,1 p,l 

W^gfWB'^, i^lbll o<^+--n/p 11/11 on/p-^ if cr + i^-n/p< 0, 

oo,1 p,1 

^II/IIb"/p-^ 115115-+- li a > -min{n/p,n/p'). 

So the result follows once noticed that B^ i ^ B^^ for any s G M. □ 
Lemma 5. Assume that a, v and p are such that 

Tl f Tl Tl\ n 

(A.22) l<P<+oo, 0<i/<- and - min - 1 < cr < v. 

p \p p' J p 

There exists a constant C depending only on p, a and n such that for aZZ A; G {1, • • • , n}, 
we have for some sequence {cj)j^z, with ||c||^i(2) = 1: 

\\dk\a,Aj]w\\LP < Ccj2~^'^\\Va\\^ri/p-i^\\w\\TDa+i^ for all ?' G Z. 

Proof. Taking advantage of the Bony decomposition (A.21), we rewrite the commutator as^ 

(A.23) dk{[a,Aj]w) = dk{[Ta, Aj]w) + dkT'^ a - dkAjTl^a . 

' ' ^ ' ^ ' V ' 

Arguing as in the proof of Lemma 6 in [12], we get 

U'-j|<4 

Now, for 1/ >0, we have 

\\VSj>_ia\\Loo < C2^"''||Va||^-. . 

DO , 1 

Therefore, for some sequence (cj)jgz in the unit sphere of ^^(Z), 
(A.24) WRjWlp < Ccj2-^''\\Va\\^-.^\\w\\^a+.. 



"^Here we use the notation T^v := T^v + R{u, v). 
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To deal with Rj, we use the fact that, owing to the locahzation properties of the Littlewood- 
Paley decomposition, we have 

Yl 9fc(5'j'+2Aj«;Aj-,a). 

Hence, using the Bernstein and Holder inequalities, 

\\Rj\\LP < C J2 II^J'+2A ^■"^11^°° ||Aj/Va||LP, 

i'>i-2 

Therefore, by virtue of convolution inequalities for series and because n/p — u > 0, 
(A.25) \\R^i\\Lp <Ccj2-^''\\Va\\^r./p-.\\w 



Next, from standard continuity results, wc know that the paraproduct and the remainder map 
B^~['^ X -Bp(^ in -Bp^^ whenever a + — n/p < and a + 1 > —mm{n/p,n/p'). We thus 
have 

(A.26) ll^fllip < Cc,-2-^'''||Va|L„/p-.||u;|L.+.. 

Putting Inequalities (A. 24), (A.25) and (A.26) together, and using classical embedding completes 
the proof of the lemma. □ 

Lemma 6. Let A(D) be a Fourier multiplier of degree 0. Then the following estimate holds 

\\[A{D),q\w\\^,T+i < C\\q\\^i-^+n/p\\w\\A.T+. 

whenever 

. (n n\ n 

u > and — mm — ,— : — 1<0"< u. 

\P P ) V 

Proof. Taking advantage once again of Bony's decomposition, we decompose the commutator 
into 

(A.27) [A{D), q]w = [A{D), T,]w + A{D)T'^q - T'^^^^^q. 

According to Lemma 2.99 in [1], we have for > 0, 

\\[A{D),T,]w\\^.+, < C||Vg||^-. \\w\\j;,.-,.. 

p,l oo,l P,l 

Next, given that A{D) is a homogeneous multiplier of degree 0, it maps any homogeneous 
Besov space in itself. Therefore the last two terms of (A.27) may be just bounded according to 
standard continuity results for the paraproduct and remainder operators. □ 
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